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SYMBOLS 


A Acceleration 

H Angular momentum of rotor blade 

w Rotational rate 

1 Moment of inertia 

T Coordinate transformation matrix 

M External moment 

ft Rotor rate 

f Rotor azimuth 

8 Flapping angle 

6 Lead-lag angle 

p x-axis component of rotation rate 

q y-axis component of rotation rate 

r z-axis component of rotation rate 

i Unit vector along x-axis 

j Unit vector along y-axis 

k Unit vector along z-axis 

e Hinge offset from rotor shaft 

m Mass of rotor blade 

R Distance to blade center of mass from hinge point 


Subscripts and Superscripts 


Fixed shaft-frame 
Rotating shaft-frame 
Blade-frame 
Principal axes frame 
Blade hinge point 
Blade center of mass 
x-axis component 
y-axis component 
z-axis component 


NOTATION 


The subscripts, superscripts and symbols used in this paper are de- 
fined below. The notation conventions are as follows: 

1. All vectors are denoted by arrows over a symbol. A symbol 
without an arrow represents a scalar magnitude. For example, w 
represents an angular rotation rate vector while w represents the 
magnitude of w. 

2. The reference point of the vector is given by the super- 

♦ *-g * _ g 

scripts. For example, the symbol w represents the rotational rate 
vector of the s' -frame with respect to the s-frame. The second super- 
script is dropped if the reference point is inertial. Instead of writing 
the angular rate of the s'-frame relative to the inertial frame as w 4 S 

. . , . -*S f 

it is then written as w 

3. The second subscript denotes the reference frame into which a 

vector is resolved and the first subscript gives the axis of that frame. 

s * 

The symbol w xs then gives the magnitude of the x component of the iner- 
tial rotation rate of the s'-frame resolved in s-frame coordinates. For 
conciseness, when a vector is resolved into the same coordinates frame 

that it references the superscript is dropped. For example, the symbol 

s * 

w xs » refers to the magnitude of the x component of the inertial rotation 

rate of the s'-frame resolved in s'-frame coordinates, so it is written 

instead as w . . 

xs 

4. The symbols p, q and r are used in place of w , w and w 

x’ y z 

respectively to represent magnitudes of the components of rotational 

s’ 

rates. The symbol w , representing the magnitude of the x component 

of the inertial rotation rate of the s'-frame resolved in s-frame coor- 

s * 

dinates, is then written as p . 

s 
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SUMMARY 


The inertial dynamics of a fully articulated stiff rotor 
blade are derived here with an emphasis on obtaining equations 
that facilitate an organised programming approach for simulation 
applications. The model for the derivation includes hinge offset 
and six degrees of freedom for the rotor shaft. Results are com- 
pared with the flapping and lead-lag equations currently used in 

the Rotor Systems Research Aircraft (RSRA) simulation model and 
differences are analyzed. 


INTRODUCTION 


A general purpose rotor model is a necessary starting point in the 
analysis and simulation of rotorcraft. In addition to allowing for 
flexibility in simulation such a model also serves as a consistent 
baseline from which a variety of more specialized models can be derived 
to satisfy specific analytical and simulation requirements. 

This paper derives the inertial dynamics of a fully articulated, 
stiff rotor blade with offset hinges and si* degrees of freedom of the 
shaft. Hingeless rotors can be represented by using appropriate values 
of spring stiffness and hinge offset in this fully articulated model. 

A configuration of current interest is the Sikorsky S-61 rotor. 

This rotor is being utilized on the Rotor Systems Research Aircraft 
(RSRA) and mathematical models are given in Ref. 1 and Ref. 2. The 
current simulation utilizes the model of Ref. 1, however the complexity 
of the model has thus far prohibited real-time simulation of acceptable 
quality on existing computational facilities. Implementation is further 

complicated by the faCt thaC Lhese equations are not presented in a 

modular, building block form that permits efficient programming. At- 
tempts to modify these equations to provide a real-time simulation on 
available computers have been hampered by the fact that their derivation 
is not documented. A lack of symmetry in some of the terms and an 
uncertainty as to the assumptions employed prompted the author oi Ref. 2 
to perform an independent derivation from basic principles as a basis 
for the analysis of computation requirements. His approach, however, 

W3t si g n ificantly different from that taken in Ref. 1, 
with those equations was possible. 


=? 


so no comparison 


The derivation performed in this paper is intended to provide a 
comparison with the rotor inertial dynamics of Ref. 1 as well as to pre- 
sent these equations in a format more suitable for structured program- 
ming. It is hoped that the results will provide a basis for more special- 
ized formulations as well as permit a more efficient programming approach 
for the general formulation. 

In Chapter 1 the derivation is outlined and the most significant 
steps are presented along with the final results. The equations are 
first derived in the blade-frame since the moments of inertia are assumed 
constant in this frame (see Appendix C). The blade-frame variables are 
then expressed in terms of rotating shaft-frame components to obtain 
differential equations for flapping and lead-lag motions. Finally, the 
rotating shaft frame variables are given in terms of the fixed shaft- 
frame rates and accelerations. Chapter 2 presents a comparison of the 
results with the equations of Ref. 1 and Ref. 2 and Chapter 3 gives the 
conclusions and recommendations. A description of the coordinate systems 
and supporting derivations are given in the Appendices. 


3 


Chapter I 

Derivation of Rotor Dynamics 


The inertial forces present in the flapping and lagging equations 
derived based on the following assumptions: 

1. The rotor shaft has six degrees of freedom. 

2. The rotor blade is hinged in both the flapping and lead-lag 
axes. 

3. The flapping and lead-lag hinges are co-located at an 
offset, e, from the shaft axis. 


4. The rotor blade is stiff. 

The equation for the rotational dynamics of a fully articulated rotor blade 
with Co-- located offset hinges, as shown in Fig. A.l, is (Ref. 3): 

1-1) M = mR x A h + H h 

The vector M h represents the total external moments 
acting about the hinges and may be written in the blade coordinate 
systems as: 


02, ^ + ^ + ^ < b 

The vector R represents the position of the blade center of mass 
relative to the hinges and is written as: 

1.3) R = R j b 

The inertial acceleration of the hinges, A h ,is given in the blade 
coordinate system as: 

K4) 

The angular momentum vector of the blade about the hinges is given in 
the blade frame as: 

1.5) 

We now substitute 1.2), 1.3), 1.4) and 1.5) into 1.1) and 

perform the indicated cross-products to obtain: 




1.6) ^ = m R A h +H h 
xb zb xb 
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1.7) 


h *h 

M = 
yb yb 


1 . 8 ) 


M h = - mR A h + H h 
zb xb zb 


Equations 1.6) and 1.8) are the fundamental relationships from 
which blade flapping and lead-lag dynamics may be derived. 

The components of hinge acceleration in the blade-frame are 
derived in Appendix D and given in terms of rotating shaft frame rates 
and accelerations as: 

1.9) 


A" = cos<$ ! . - e(r , - q ,p ,)| - sin & , A* , - eCp^, + 

xb - xs T s’ s 1 s’J : ys' K s' M s J 


2 O 


1 . 10 ) 


A yb “ sin6 ' A zs' + e(p s' + W 5 ' 


+ cosb sin<5 A® gf - e(r gI - q s ,P g ,)| + cos6 [ A ys' " e( P s * + r s P 


2 >] 


1 . 11 ) 


h 


A = cosb'.A , + e(p , + q ,r 

2b *7C f r c. 1 c ’ 


s’ s 


r 


+ sinMsin* 5 A S , - e(r , - q ,P ,) + cos^iA S , - e(p 2 , + r ,), 
| xs' s s t- ys' r s' s 

The blade-frame components of the time rate of change of the 
angular momentum of the blade about the hinge are derived in Appendix C and 
are given in terms of rotating shaft-frame rates as: 

1.12) H' 

+ 

- sinB cosB- 


1.13) 


1.14) 


• , p * • 

xb = V.’ 008 * " q s ' sin6 “ < p s’ Sin6 + q s' f os6 > 6 * 3 
(cos 2 3 - sin 2 8)j(p gl sin<5 + q g ,cosS)(r g , - ^ 


1 |cr s f - b 2 - (p g1 sin<5 + q gf cos<5 ) 2 ^ | 


.h 


H y b “ I yy\ ■Ps * sin6 + q s « COS(S + (Pg.costS - q g i sin6) 5- (r g ,-6)g cos 


I 


r ,-b + (p ,sin6 + q ,cos£)3 sin? 
s s s 


[. 


H* 1 = I, sinb-jp , sin<S + q^, cos<$ - 2(r , - b) 8 + r , (p ,cosb 
zd d i ^ s ^ s as 

- q . sin<$); + cos8', 26(p ,sin<$ + q cos^) + r . - & 

s \ [ s ‘s' s 

- (p T sin5 + cj , eos^)(p t cos6 - q jSirnS) 1 ' 

s s s s \ J 
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Substituting 1.11) and 1 . 12 ) into 1.6) and solving for g gives: 

l - W ‘ ‘ V + + VV>) 

+ slnB [si n s jA M , - e(r s , - 9 S ’P S '* + cosS jd ys t 

- e(p s , 2 + r s 2 )]| 

+ P s *cos<5 - q s ,sin6 - (posing +q ,cos 5 )6 


S -i S - '«- s n s 

+ (cos 8 - sin 2 g) [(p s ,sin 6 + q s ,cos<s)(r , -$)> 


| vr 's' ' I s i x * s 1 

sing cosg . (r g , - 5 ) - (p s ,sin <5 + q ? ,cos^) 2 » - M x h 


Substituting 1.9) and 1.14) into 1 . 8 ) and solving for 6 gives 


1.16) 


6 " I b cos 8 | cos6 L A xs' - e <V ' VPs' 5 ! ~ si "6 
2 ' " “• * 

3 

s 


r 


-e(pf, + 




+ ta ” ^s’ 31 ” 6 + V co3 « - 2<V -4)5 + r s ,(p s ,co 8 4 
-q s .si04)l 

+ ^B (p g iSin4 + V <=° 3 » + f 3 , - (p g , si„6 + % , coss) (p ^ coa( 


- zb 


- q t sing ) 


I^cosg 


H.e rotating shaft-frame components of rates and accelerations are 

derived in terms of fixed shaft-frame components in Appendix B. The resulting 
equations are: 


1.17) 


p s . = P s s ii# + q g cosip 
V = cos ^ + q s sin ip 


V " r s~ Q 


1.18) 


P s * PgSini/) + q cosi|/ - q ,!J 


q s’ = * PgCosi*, + q s sint|f + P q 

• • 

V = r s - 
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1.19) 


A ,= A sin'*' + A cos* 
xs xs ys 


A , = - A cos* 4- A sin* 
ys xs ys 


A , = A 
zs zs 


Equations 1.15) through 1.19) then allow determination of flapping 
acceleration and lead-lag acceleration from fixed shaft-frame components 
of rates and accelerations. 
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Chapter II 


Co mparison with RSRA Rotor Equations 


In this chapter the results of Chapter I are compared with the equations 
of Ref. 1 and Ref, 2 to evaluate the relative merits of the equations for 
simulation and analytical studies. The blade flapping and lead-lag equations, 
as given in Ref* 1, are: 

2.1 ® cosBa 4 e 2^(P cosi* - q sino) 4 p sinC + q cost' 

I K z s r s y M s Y J-s M s y _. 

2 


4 sinScos ' A sinijj - A cosil; - e(r - ft) 
ys xs s 

2 * * * 

4 cos 6 cos6 p^sin 1 ^ 4 q^cost^ - 2(6 4ft) (q^sini^ - p^cosi^) 


- 

2 xb 


2.2 


- 2ftsin6 p sin^ 4 q costJ> 

t S b 

4 cosgsinB 26 (r - ft) - (r - ft) 

s s x, 

— D 

TO p O 

* *; r sin6 A sint - A cosi l> - e(r - S2) ' 

I-cosB i ys xs s l 

o r i ^ # “1 i * 

- cosSA sint + A cost + e(sl - r )r 

xs ys s | * 

+ tan 6 28(ft 4 6 — r ) 4 q sin (’4/ + 6) - p cos (\p 4 6) 

^ S S S „ , - 

+ r - ft4 2gjcos6(q sinft - p cosft) 4 sin6{p sin*J> 4 q costy) 
s ^ s s s s 

< 

- zb 


I cos 3 
b 


where the notation has been changed as required for consistency with this pape 
In order to compare equations we must write g and 5 in terms of s-frame 
variables. This is accomplished by substituting 1.17) through 1.19) into 
1.15) and 1.16) to get: 

** mR I ^ 

2.3) 8 * j — !cosg A 4 e 2ft(p cos\p - q sin^*) 

i, t _ z s s s 

b . 

, * 

m * 

+ p sinip 4 q cosij) 4- r (q sin^ - p costi)' 
s _ s s s s 

r 7 2 

+ sin£cos6 A sint ~ A cost - e (r - i]) _ + (p sinf + q cost) 
ys xs L s s s 

4 sin6sin6 i A sinijj 4 A cosip - e :(r - ft) 

xs ys u ^ 

- (p sin^ 4 q cosj)) (q sim|> - p cos\j>) ^ 

S S-S S _ 1 1 

2 2 ♦ ~ 

+ (cos 8 “ sin 8) cos6(r - ft -6)(q sini* - p cos^) 

, • ' s s 


. 1 \ 
4 sin6(r - ft - 6) (p sinip 4 q cos^)- 
s s s 


+ singcosB . (p sint + q cest)sin6+(q sint-p cost)c os 3 -( r - ft-f) 
I-.S s s s *^s 

4 * • 

sint + q cost - 0* + £)(q sint - p cos r ) 

s S S S i 

- sin 5^ 


sin^ 


p costJj 4 (ft 4 6)(p sin'^ 4 q cos^ 4 - y.b 

s s \s — 

* x b 


2,4) 6 I^cos 8 [ sinS ^ys Sin A xs cos ♦ 

~ 2 j 

- e (r g -ft) + (p^sintjj + q g cos^) 2 J 

- cos<5<A sin^ + A cosii.' 

xs ys 

r • 

+ e [fl ~ r g + (p g sin^ + q g cos<j,) (q g sin,J, - p cos^)j <■ 

r . . s J - 

+ tan6 | 2&(tt + 6- r g ) + q g (cos^sin6 + sin^coso) 

- p s (cos^cos<5 - sin^sin5 ) 

r s {^Pg (sinJp cos<5 + costpsind) + q g (cos^cos6 - sin^sin6)]j 
+ r s " ° + 26 [cos6 (q s sin^ - p s casp) + siivS Cp g sii^+ q g cos^)] 
+|^ s (co^cos6 - sinpsind ) - q^cos^sinS + sinj,cos6)J [_p s (sim/,cos< 

+ cosifisinS) + q g (cos^cos6 - sin^sin6)J 

M h . 

- zb 


I. cos 6 
b 


Neglecting second order terms in shaft rotational rates (p , q , r ) and 

s n s s 

using the trigonometric identities: 

2.5) sin(ty + 6 ) = cousin 6 + sin^cosd 
cos +6) = cos^cosg - sin^sin <5 
we may rewrite 2.3) and 2.4) as: 


2 . 6 ) 


6 “ I h [ co ^J A zs + e ^2^(P g cos4) * q g sin^) + j^sin^ + q^os^jj 

+ sin3cri«6.C\ c _ a -x 2l 


+ sinBcosfijA^sin* - A^cos* - e(r g - a ) 
+ sinBsin5 f A xg sinij; + A yg cos^ - e(r’ o - ^ 


2 <■“ r ys r ' 8 ij 

+ (cos s - sin 6) cos6[- ( n + 6)(q g sinif - p cosj.)] 
+ sinS^- (ft -f 6)(p g sin4* + q cos^)fj 


G) ~ (r - ft)' 
s 


--) 


+ r . 

cosgsing ( 26(r - 

r |_ S V‘ s «/ - o r 

+ cos« lyin* + q gCOB * - (ft + ^[q^in* - p & cos^ | 

- sinfip s sinp - PgCoaji + (ft + 6)[p gS in^ + q^os^jj- 


M 
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2.7) 6 = ^- R 


2 , 


I b costi |_ sln5 A ys sln ^ - A xs cos ^ - e (r g - fl) 

r 'O 

- . . . . 1 1 

- cosi A sin{i + A cos* + e(i> - r ) 

I aS yS g 

r . . . ^ 

+ tanji » 2{5(id4cj - r ) 4 q sin(p f 5) - p cos(u 4- 6)> 

L & S s i 

• • » 

+ r g - ':+ 23ICOS6 (q g sin* - p cos* ) + sin5(p sin*+ q cosij;) 

M h S 3 

- M zb 


I.cosg 

b 


Comparing 2.7) and 2.2) we see that the equations are identical so the 
only additional assumption made in the lead-lag equation of Ref. 1 is that 
second order terms in shaft rates may be neglected. 

Comparing 2.6) and 2.1) we see that further rearranging is required for 
comparison. Using the identity: 

2 . 8 ) 


2 2 
cos 8 + sin 8=1 


we multiply the last two lines of 2.6) by 2.8) and regroup to get: 


r 


2 ’ 9) 6 " I I COs6 i A zs + e ! 2ft (p cos* - q sin*) + p sin* + q cos*j 

b b s s -J 


+ sin8cos<5<A sin* - A costf/ - e(r - o) 2 ' 

J ys xs s 1 

A simp + A cos* - e(r - * vJ 
XS yS r s 


2 F | 

+ cos e. lC os6|_p s siiu|» + q cos* - 2 (fl + 6)(q sin* - p cos*)] 
I s s — I 


l 


L 


*)!' 


sin^L^gSin^ P g Cos^ 4 2(Q 4* 5)(p sin^ 4 q cos 

r. s i s i 

+ cosSsing <26 (r - ft) - (r - ft) 2 - 6 21 
2 p ^ s s J 

+ Sln 2 e f cosS^sin* + q^cos) - sin«[q sin* - p cos(j> 

H h L J 

- X b 
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A comparison of 2.1) and 2.9) now reveals the following discre- 
pancies in 2.1) 

1. Line 3 in 2.9) has been neglected. This may be justified 
by assuming sin 6 is small. 

2. All terms not including the factor fi (rotor rate) have 
been dropped from line 5 in 2.9). This may be justified since the other 
factors are relatively small compared to A but the same is true for line 

4 and similar terms have not been neglected there. 

->. The square of the lead-lag rate, 6 , has been neglected in 
line 6. It is not clear however, that this is negligible compared to squares of r 


and products of r and <5. 

s 


4. Line 7 in 2.9) has been neglected. This may be justified 
? 


by assuming that sin^B is small. 

The approach of Ref. 2 is to substitute Euler’s equations in the blade- 
frame into the moment equations and solve for the rotational accelerations 
xn the blade-frame. In this paper, Euler's equations in the blade-frame are 
given in equation (C.3) and the moment equations are equations 1.6) through 


1.8). Performing the substitutions we obtain: 


2 . 10 ) 


K b I, 


<M xb ~ mRA .K> " 9u r 


zb' 


“Vb 


2 . 11 ) 


r . — r (M . + mRA* 1 , ) + p q 
b I fe zb xb H b 4 b 


Equation 2.10) agrees with its counterpart in Ref. 2, however the equa- 
tion for r b in Ref. 2 has a minus sign on the hinge acceleration term and 
appears to be in error. 
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In Ref. 2, p fa and r fe are integrated and the flapping and lead-lag 
rates obtained as a function of the difference in blade-frame and rota- 
ting shaft-frame rates. The appropriate equations in this paper are given 
in (B.16) and the resulting expressions for flapping and lead-lag rates are: 

2.12) = p ,cos6 - q f sin5 - p 

s s b 

2.13) 5 = [p sind + q .cosdlsing + r ,cosfc> - r 

L- S S i s b 

No differential equation for has been obtained since the moment 
of inertia about the y-axis has been assumed negligible so there is a 
problem in obtaining values of for use in 2.10) and 2.11). In Ref. 2 it 
is suggested that, with the exception of blade pitch control, may be 
negligible. This seems unlikely, however, since the rotational degrees of 
freedom of the rotor shaft as well as the lead-lag rate (5) and the rotational 
rate (12) can couple into the y-axis of the blade-frame. A solution to this 
problem can be obtained by noting that 2.12) and 2.13) have been obtained 
from the p^ and r^ equations of (B.16) so an independent equation for 
remains and is given as: 

2.14) q^ = Jp^.sind + , cosS j cosh - (r g , - d) s i n g 

where 6 is obtained from 2.13). The hinge accelerations in the blade-frame 
and the rotating shaft frame rates are required in 2.10) through 2.14) and 
may be obtained from the fixed shaft-frame rates and accelerations by 
1.9) through 1.11) and 1.17) through 1.19). 


Chapter III 
Conclusions 


The equations for the inertial dynamics of a rotor blade as derived 
in Chapter I are intended to serve as both a general purpose simulation 
model and as a baseline for more specialized analytical studies. In this pa- 
per, no assumptions as to the relative magnitude of the rates have been 
made. For most applications it should be possible to further simplify 
these equations by appropriate assumptions. This derivation has been in- 
tended to serve as a baseline for more specialized applications; hence, only 
the most generally applicable assumptions have been utilized. 

By writing the flapping and lead-lag equations in terms of rotating shaft 
frame rates and accelerations the equations have been significantly simplified 
over the formualtion in Ref. 1 where these equations were given in terms of 
fixed shaft-frame variables. In this paper the rotating shaft frame variables 
are obtained by a separate transformation from the fixed shaft frame variables 
This allows for a more structured approach to programming and thereby improves 
the efficiency of the simulation. Programming errors are also more readily 
located and corrected in this format. 

The validity of the equations of Ref. 1, which are used in the current 
RSRA simulation, was checked by combining the modules derived in Chapter I 
to obtain differential equations for blade flapping and lead-lag motion as a 
function of fixed shaft frame variables. By neglecting higher order terms in 
some rates the lead-lag equation was f , ^ .o De identical to the lead-lag 
equation of Ref. 1, but discrepancies were found in the flapping equation. 

The discrepancies appear to be the result of inconsistencies in applying small 
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an 8 l. and low race approximations in the equations „ f ^ ^ 

haa been made an this paper to determine the validity of these assnmptions 

or the numerical significance of the discrepancies. 

The approach used in Ref '? „ K r • 

er. . is to Obtain equations for the rotational 

accelerations of the blade-frame and integrate these 

■ Lnce g rate these accelerations. The 

flapping and lead-lae rate<; av- D eh <- 

la* are then found as a function of the difference 

Between blade-frame and rotating shaft-frame rates. Phis approach was re- 
peated in Chapter II for comparison purposes and a discrepancy was found in 

eq “ atl °- *■ Si8 " ° f — acceleration ter. in this equation 
appears to be in error in Eef. 2. 

blem with the approach taken in Ref. 2 i s that no differential eq,« 

or q b can be obtained since the y-axis moment of inertia of the blade 

a assumed negligible end the suggestion made in Kef. 3 that „ b may be email 

and can be neglected does not appear to be justified 

J tlIled - A solution to this 

problem is suggested in Chapter II of 

P °f thlS pa P er and involves solving for q 

function of rotating shaft axis rates and lead-lag rates. With this 

addition, the equations of Kef. 3 M y be well sulted „ ^ 

BBey are Presented in an extremely moduiar format. Por analytical worh 
Bowever. explicit differential equations for blade flapping a „d lead-lag 

motions are required rather than expressions for these verfables in terms of 

the output of other differential equations Th 

he equations derived in this 

Paper are intended to satisfv , , 

satisfy both simulation and analytical requirements. 
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Appendix A 
Coordinate Systems 


The four coordinate frames used in this analysis are shovn i„ 

A. 1. They are: 

1) The fixed shaft-frame (s). 

2) The rotating shaft-frame (s'). 

3) The blade-frame (b). 

4) The principal axes frame (p). 

The fined sheft-fr* (s, is centered i. the rotor huh eith its 

X ^ ^ ~ Pl “' -* «. 2-axis normal to the 

plane of the rotor hub. The s -fr». r 

S th * startin 8 Point of this analysis 

Rotational rates and accelerations and translation*! 

a translational accelerations of this 

frame are assumed available In terms of rates an H 

rates and accelerations at the 

aircraft center of gravity. 

The rotating shaft-frame <»•) is shoes relative to the fixed araft- 

1,1 Fl8 - A ' 2 - 1818 f -“ ™ in the huh but rotates with 

the blade. Its y-axis Is directed through the co-located hingea and Its 

”~ 1 to the rotor hub - - — - - — - — 
with the y-axis of the s-frame at * = o. 

The relative orientations of the blade-frame and the rotating shaft-fra 

are shown in Fie a "! m, 

8- A. 3. The rotating shaft-frame is rented . . 

hinges and fixed to the blade " 

, . ’ S ° tS ° rl8in 15 ^.placed from the s-frame 

origin by the amount of hinge offset. The y-axis of th. r 

with the blade and its x-axis is , **** 

normai to the plane of the flapping hinge. 
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The origin of the principal axes frame (p) is co-located with the 
origin of the blade-frame but the principal axes are fixed in the blade 
and oriented to eliminate all products of inertia. 

The relative orientation of the principal axes frame and the blade- 

frame is shown in Fig. A. 4. It has been assumed that the y b axis (axis 

of feathering) passes through the center of gravity of each blade cross 

section. The y b axis is then a principal axis since the x-y and z-y 

products of inertia vanish. The x and z principal axes are consequently 

m the x-z plane of the blade-frame at an angle 9 from the v and z 

P D b 

axes, where 0 p is chosen to make the x-z product of inertia vanish. 

Note that 0 p is a function of the blade twist and the impressed blade 
pitch, so it will vary with cyclic and collective control inputs. 

The transformation from the blade-frame to the principal axis-frame 
consists of a single rotation, 0 p , about the y fc axis and is given by 

~i 

cos0 o -sin0 

(A. 1) T , = 

p/b 0 1 0 

sin© 0 cosQ 
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Appendix B 

Rotational Rate of B1 ade Frame 


Sorting from the fixed shaft— framo u 

shaft frame, s, the rotational rate of the 

rotating shaft frame, s ’ , is given by: 


(B.l) 


±s 


S s ^ ^s'-s 


W = w a + W 


where w s is the rate of the s frame ^ f 

given in s-frame components a 


(B.2) 
The vector w 


w 


;s'-s 


P s^s + q «X + r it 
s s S s s s 


and from Fig. A . 2 is seen to be: 
ts'-s 


1 S the rate of the s '-frame relative to the s -frame 


(B.3) 


■-0 k 


where fi is the rotor rate. 

Combining (b.l,, (B.2) endfB.B, the rotational rate of the o' frame 
nay be written as: 


(B.4) 


s' 

w * p 


“ S ' 8 ' t_ + q 8 ’ * ^ -s' * 


S S 


J + r 


s \ 


where : 


(B.5) 


s f 

P = P 
s P s 

s' 


r = 


r s ^ 


From Fig. A. 2 the transformation from the shaft to th« 

ne snar t to the rotating sha 

axes is seen to be: 


(B 6) r 8lrf co ^ 0 

^ T sVs " ‘ cp ^ siltf 0 

0 0 1 

Me 

where ip is the azimuth of the blade. 


The rotational rate of the s' -frame may then be resolved into s* 
components using (B.6). We then have: 


(B.7) 
where : 


-V e * 

v s 


V 


V + v V + v V 


(B.8) 


P g » » P s sint(» + q g cosiJ> 


P s cosij) + q g sin^ 


r • ■ r - o 
s s 


Differentiating (B.8) the rotational accelerations in the s'— frame 
are found to be: 

V- P 8 si »» + 4 S cos* - q g , n 

(B.9) q gt « -p g cos* + q g sin* + p gl Q 


V - r - a 


Using (B.6) the acceleration of the origin of the fixed shaft-frame 
may be written in rotating shaft— frame coordinates as: 

A XS' = A XS Sin * + A ys COS « 

£ 

(B.10) ky s r -A xs cos\fi + A ys sin* 

A S s 

A , = A 
zs zs 

The rotational rate of the blade— frame , b, is given by: 


(B.ll) 


w* 3 , ws' + wb-s* 


where w is the rate of the blade-frame relative to the rotating shaft- 
frame* From Fig. A. 3 we may write w^ s in blade coordinates as: 

(B.12) ^' S » - 6^ + 6sin6j b - 5cos8 ^ 

where 8 is the flapping angle and 6is the lead-lag angle. 
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The transformation from rotating shaft-frame to blade-frame is 
found from Fig. A. 3 to be: 

— , r n 

1 0 0 j ] coss - sins 0 

0 cosg -sing !j sinS coss 0 

0 sing cosg j ■ 0 0 1 


or 

(B.14) _ 

Vs’ = 

Usin" (B.14) to resolve (B.7) into the blade-frame and combining 
the result with (B.ll) and (B.12) vre get: 

(B.15) « b * P b h + % \ + r b \ 
where: 

p b * P s t c °s6 - q gl sin« - 3 

(B.16) * jjp^.sinS + q g , cosfiJcosB - (r g , -6) sinB 

r b * [fs* sin6 + q s' cos< D sine + ( r g i “h cosB 


cos 6 -sinS 0 

cosBsinS cosgcosB -sinB 
sinBsinB sinBcosB cosB 


(B.13) 


l b/s' 


P b = P s « cosS _ <? s . sinS -3 - Ip, sinS + q g , cosSJ 5 

(B.17) q, = ip .sinS + q , cosS + Ip C os6 _ q f sinfi | g 
- (r , - 6)&lcosg 


f. " 

~ , - 6 + p . 


161 A sii 


, g . U> s »sin5 + q g , cosSj glsing 

• f • ^ • i— J — 

r b * (P s » sinS + q g , cosS + |p gf cosS - q gl sinSj 5 


‘ ( V 


-6)3 • 

*• - 


sing 


r* 

+ r s' " 6 + [p s » sln5 + <l 8 i cosSjglcosg 

m* 
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Appendix C 

Euler's Equations in the Blade Frame 


The angular momentum vector in the blade— frame may be written in terms 
of the principal moments of inertia as (Ref. 3): 


-»-h T 

"p/b Vp/b W b 


(c.i) = r „ i t „ w, 


where T ^ is the transformation from the blade-frame to the principal axis 
frame, as given by (A.l), 1^ is the principal moment of inertia tensor, and 
^ is the angular velocity of the blade-frame, as given by (B.15). 


Substituting (A.l) and (B.15) into (C.I) gives: 




H 


'xb 


2 2 

(*xx COS + s '*‘ n Pi. + " I„„) sin0_cos0_r. 


zz 


r b 


zz xx 


Pb 


(C.2) 


ljb r 

H , = I q 
yb yyb 

H *h = (1 2Z ' I xx ) sin0 « cos0 «Pk + sin 2 0„ + I.„ cos 2 0_)r, 


P r b 


XX 


zz 


P' b 


The time rate of change of the angular momentum vector is then: 

(C.3) H h = H h_b + w x H 

D 

A significant simplification can be obtained by the following assumption 

A 


(C.4) 


I = I S L 
xx zz b 


The angular momentum components are then: 

H h = Ip 
xb b K b 

(C.5) 5 h , « I q, 

yb yy M b 

* I r. 
zb b b 

Note that this assumption has eliminated all 0^ terms, making the moments 
of inertia constant in the blade-frame. 
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Using (C.5) in (C.3) gives: 


H 


xb " Vb + (I b" ^Vb 


(C.6) 


H 


I q, 

yy H b 


zb "b'b 

If we further assume: 


yb 

s zb ■ Vb + -y Vb 


(C.7) 


I << I 

yy bb 


then (C.6) becomes: 


o , 

xb 


(C.8) 


V p b + Vb } 


H n - T • 

H yb - Vb 


*b " I b (r b- p b q b ) 


Using (B.16) and (B.17) in (C.8) gives: 

“xb " V, V cos6 ”*i s » sl “«S- [p s t sinfi + q^coss] $ _e 
[c°s 6 ~ sin ^[p s * sin5 + q gf cosfij (r.,-6 ) 

f < t. 2 _'l 


(C.9) 


1 JL a *s “J ' s»' u ' 

-sine cose [(r s ,.J ) 2 -<p s , slnS _ cos{) : 

r 


(C.10) 


(C.ll) 


yy 


^Pg.sinfi + q s ,cos6 + (p s ,cos6- q s? sin6)5-(r s , -6 
’[ r s* " 6 + (P s -sin6 + q s ,cos6)^ sing) 


H h = T 

zb “ *b 


+ cosg 


jsing [ p o , 

H 


J 


oin6 + 5 s ,cos6 -2(r s , -j)e + r>t f p ^,co: 

4 -_ *1 , * *’ L 

J + V - 


L r s 

Pg.sinfi + q gf cos6 


L - i gf - c 

-[p s .sin« + q s ,cos«j f Ps , c „s« -, s ,si„6]l 


cosB 


- 

q s ,sin6 
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Appendix D 

Acceleration of Hinge Point 


The acceleration of the hinge may be written as: 

(D.l) £ h = I s + t h ~ S 


where is the inertial acceleration of the origin of the s-frame 
and 1 S is the acceleration of the hinge relative to the s-frame. 

Using the transformation of(B.6) we may write interms of its 
components in the s '-frame as: 

(D.2) V 


ts a s T , ,s -> . . s -> 

A xs’ i s’ +A ys’ j s + A zs' V 


From Fig. A.l the acceleration of the hinge relative to the s-frame 
is seen to be: 

(D.3) A* 1-3 - i 

where: 

(D.4) e »e] 

o * 


We may write the derivative of (D.4) in terms of the rotational rate of 
the s 9 — f rame as : 

(D.5) e * e *5 gf + w S x e = w S x 

where we have noted that e is zero in the s’ -frame. The second derivative 
is then: 

(D.6) "e * w S x "e + x e* 
combining (D.l), (D.3) and(D.6) then gives: 

(d.7) l h * t s + i s ' x % + n*' x x %) 

Using the definition of w s in (B.7) and its derivative along with 
(D2) and (D.4) we obtain: 


(D.8) 
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where : 


A , = A s , - e(r r - q , p T ) 
xs xs f s M s s 

(D.9) A h . = A S , - e(p } + r , 2 ) 
ys ? ys r s* s 

A zs' = A zs’ + e( P S ' + V V* 


Using the transformation of (B.14) we may then write A n in terms of 
its components in the blade-frame as: 


(D.10) 

Ah = 4 S + A Hb \ + A zb \ 


where : 



D.ll) 

A*) * A* 1 »cos6 - A* 1 t sin6 

xb xs ys 

A yb * (A xs • sln6 + A y S * cos6) cos B - A^ s , 

sing 


A zb " (A x a f Sln<5 + A J S » cos< 9 sin6 + A^, 

cos6 
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